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We report on a lattice calculation of the form factors Aq and A3 for the pseudoscalar 
to vector meson semileptonic decay p'^l~v'i. We find that resonant (or pole-type) 

Oh! contributions alone are unable to describe these two form factors simultaneously. For 

^ I the quantity AQ{q'^=0), which is important phenomenologically for the determination 

of \Vub\, we extract a range of values, Ao((?^=0) = (0.16-0.35)lg, where the range is due 
to systematic uncertainty and the quoted error is statistical. We have also determined 
>< ; A,iq'=0) = 0.28tltt. 



Abstract 



Introduction 

The semileptonic decay p^l^vi is determined by the matrix element of the V—A weak 

current between a B meson and a p meson. The matrix element is, 

{p{k,r^)\u^^{l-^,)h\B{p))=^*PT,p, (1) 

with form factor decomposition, 

T^p = ^^^^ ^ ef^^spp^k^ - i{mB + mp)Ai{q'^)g^p 
niB + frip 



niB + mp 



where q = p — k is the four-momentum transfer and rj is the p polarisation vector. The form 
factor A can be written as 

Aiq')=Aoiq')-Asiq'), (3) 

where, 

Znip ZTUp 

with Ao(0) = ^3(0). In the limit of zero lepton masses, the term proportional to A in 
equation (0) does not contribute to the total amplitude and hence to the decay rates. Pole 
dominance models suggest that V, Ai for i = 1,2,3 and Aq correspond to 1~, l"*" and 0~ 
exchanges respectively in the t-channel |T]]. 

Neglecting corrections suppressed by inverse powers of the heavy quark mass M, the 
following relations hold when is close to the maximum recoil value g^g^,^ = (M — m^)^ 

= const, AsG/M^/^ = const, (5) 

where arises from perturbative corrections and is chosen to be 1 at the B mass. In leading 
order 0, 

e.e(M.™«)4i^)*. (6) 

In the calculations reported below, we will use /5o = 11 in the quenched approximation and 
Aqcd = 200 MeV. 

The differential decay rate dT{B'^ — > p^l^ui) / dq^ at = is given byQ, 



dV{B^ ^ p+l-ui) Gl\Vui,\'' 



dq^ 



q2=o 19271^171% 



{ml-mirAo{0)\ (7) 



and is determined by the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix element \Vub\ 
and the form factor y4o(g^=0), which accounts for the nonperturbative QCD contributions. 
A lattice measurement of Aq{0) can therefore be combined with experimental measurements 
of the differential decay rate to extract \Vub\- 

^The decay rate expression in equation is naturally expressed in terms of |^3(0)|, but since our lattice 
data for Aq are of much higher quality, we prefer to use the relation Ao{0) = ^3(0) and express the result 
in terms of |ylo(0)| 
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It has previously been suggested || that hadronic uncertainties can be reduced by 
combining a measurement of the differential rate with the recent measurement of the rare 
radiative decay B K*'y to determine \Vub\, according to the relation: 



R{B ^ K*-f) 



B 



m 



K* 



F iKfeP 



[m 



B 



m 



2A3 



mi 



Q^r-pB^K* 



(0) 



,B^p 
^0 



(0) 



Here, the hadronisation ratio R{B — > K*^) = r{B K*'j)/r{b 57) depends only on the 
form factors 71(0) = iT2(0) for the exclusive radiative decay, and is given up to (9(l/mf) 
and perturbative QCD corrections by,0 



R{B K*-f) 



\miy J 



m 



ml 



7^2(0) 



(9) 



■^B^K 



(o)Mr 



^''(0) is 
Lattice 



Using equation (H) to determine \Vuh\ is advantageous if the ratio 2T2 
known in advance. There are some predictions that this ratio is close to 1 ^ 
studies can help here by allowing a model-independent calculation. 

QCD sum rule analyses - |T5[, quark models []16|, ^\ and a combined theoretical and 
phenomeno logical analysis [|18| give the following picture for the form factors for B ^ p 
transitions: V and Aq have significant dependence, while Ai is rather flat in and A2 
shows intermediate behaviour. 

In this note, we present the first lattice study of the form factors ^0(9^) cind ^43(5^), and 
give a range of values for AQ{q^ = 0) which, given experimental measurements of dr{B^ — >■ 
p'^l~h'i)/dq'^\q2=Q, may be used to determine \Vub\, according to equation We first study 
Ao(g^) and then argue that a combined fit to Aq and A^ can give more accurate results 
for Aq{0) as lattice data improve. We also report on a value for A2{q'^=0). Other lattice 
calculations of A2 for B^ p^l'vi can be found in references PUI ETl 



Lattice Details 

The results described below come from 60 SU{?>) gauge configurations generated by the 
UKQCD collaboration on a 24^ x 48 lattice at /5 = 6.2 in the quenched approximation. The 
0{a) improved Sheikholeslami-Wohlert (SW) |^ action was used for fermions, with "ro- 
tated" fermion fields appearing in all operators used for correlation function calculations p4 



Three-point correlators of the operator in equation (|T]) with a heavy pseudoscalar meson 
(the "5") and a light vector meson were calculated, for various combinations of heavy 
(p) and light (k) meson three momenta. We refer to each combination with the notation 
|p| |k| where the momenta are in lattice units of 7r/12a (a subscript ± on k indicates 
that p and k are perpendicular). Here, the inverse lattice spacing determined from the p 
mass is = 2.7(1) GeV We use this value for since we are dealing with decays to 
a p meson, but it should be remembered that different physical observables lead to different 
values for (UKQCD determinations of are summarised in [Q). This uncertainty in 
should have little effect on the dimensionless form factors evaluated in this paper. 



^The theoretical prediction for this ratio may be subject to long-distance effects [0, || 
■^The CLEO collaboration have measured exclusive semileptonic B ^ p and B ^ tt decays, and have 
preliminary results for distributions p9|. 
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channel gVGeV^ Ao{q^) Asiq"^) Af'^iq^) 



1 - 


^V2i_ 


14.4li^ 


0.98t\l 




-0.501?^ 


- 


.V2 


15.3li 


OMtl 


-0.52i^? 


-0.64^1 


1 - 


l± 


16.7ti 


1.26^1^ 


-0.34li 


-0.25^^8^ 


- 


1 






-omttl 


-0.64^1 


1 - 





19.71^1 


2.24^^^ 







Table 1 Values of the form factors AQ{q'^) and ^3(5^) for p'^l vi. Errors are statistical 

only. The final column contains A3 values obtained from previous results for Ai and A2 [30|. 



For relating lattice results to the continuum, we have used the perturbative value for the 



renormalisation constant of the axial current operator p7 | 



0.97. 



(10) 



Uncorrelated fits were used for extrapolations in the heavy quark mass, although the ex- 
traction of the form factors from calculated three-point correlation functions used correlated 
fits |]2S|. Statistical errors are 68% confidence limits obtained from 250 bootstrap samples: 
unless otherwise specified, errors quoted below are statistical only. More details of the lattice 

^ and 
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calculations can be found in references 

The SW improved action reduces the leading discretisation errors from 0(a) in the Wilson 
fermion action to ©(a^a). For quark masses m around that of the charm quark, a^ma can 
be of order 10%. For this first study we have not tried to disentangle these m-dependent 
discretisation effects from physical mass dependence. We leave them as a component of the 
systematic error. 



Computing Ao(0) 

We have extracted values for Ao(g^) for five momentum channels (we cannot extract Aq for 
the zero recoil channel, which has q^ = Q'^ax)) f*^^ f*^^^ heavy quark masses around the 
charm mass. We choose channels where the value of lo defined by, 

, M^ + m^-q^ 
00 = V ■ V = — , (11) 

is nearly constant as the heavy quark mass varies, where M, v and m, v' are the masses 
and four-velocities of the pseudoscalar and vector meson respectively. For these channels, 
we have extrapolated linearly and quadratically in 1/M to the B scale at constant tu, by 
appealing to heavy quark symmetry. This gives five values for the desired heavy-to-light 
form factor Ao(g^) with q^ values near g^ax- they are shown in table (the table also gives 
values for A^[q^) for later use). 

These values can then be fitted to various q^ dependences for Aq and an extrapolation 
made to determine Ao(0). We will refer to this as the "constant cj" method below. As a 
check, we also extrapolate to = for each of our four heavy quark masses around the 
charm mass, and then extrapolate Ao(0) to the B scale: this will be referred to as the "g^ = 



extrapolation" method. Our methods are described in more detail in references ||2£] and . 
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no 


q^=0 extrapolation 
^(0) 


Ao(0) 


constant u 
mo / GeV 




1 


0.35t^ 


0.341^ 


4.8l? 


0.32 


2 


0.2411 


0.231^ 


5.41^ 


0.39 


3 


o.ieti 


0.191^ 


6.01^ 


0.42 


4 


0.091} 


0.171^ 


6.51^ 


0.44 


5 


0.041} 


0.161^ 


7.01^ 


0.45 



Table 2 Values of the form factor Ao(0) and the mass parameter mo for — > p^/ vi for different 
assumed dependences given in equation (p^). Errors are statistical only. 



We have fitted Ao{q^) to dependences of the type, 

Ao{q^) = "^'j^^, , (12) 

for a range of "polar powers" no = 1, . . . , 5. The higher powers, while not physically mo- 
tivated, allow us to explore different dependences away from the pole region. The 
dependence of Aq in equation (^) is consistent with the heavy quark symmetry requirement 
of equation (^, which applies for close to g^ax^ provided that. 

This dependence on the heavy mass has been used when scaling Ao{0) values from the 
charm mass to the bottom mass in the g^ = extrapolation method (see figure |I|). All 
extrapolations in the heavy pseudoscalar mass M have been performed to (9(1/M^). The 
results for both methods are given in table ^. 

Table |^ shows that the g^ = extrapolation method gives results compatible with the 
constant uj method for powers 1, 2 and 3 in the fitted g^ dependence, but that the agreement 
gets worse for higher powers. This discrepancy is not significant however, because the 1/M 
extrapolations are not under control for higher polar powers. Indeed, for the lower powers 
the linear and quadratic 1/M extrapolations agree, but for the higher powers they do not, 
indicating that higher 1/M corrections must be included. This point is illustrated in figure |1]. 
For the constant u method, in contrast, the linear and quadratic 1/M extrapolations all agree 
within errors for each momentum channel studied. Our results for fits to v4o(g^) using the 
extrapolated points at the B scale are shown in figure 

In summary, we find that there is some uncertainty in the extrapolation to g^ = (in 
contrast to our finding for Ai in our earlier studies |3^). The uncertainty is smaller for the 



constant uo procedure, but nonetheless, a broad range of 0.16-0.34 has to be given for Ao{0). 
Higher powers of no do not increase this range when Ao{0) is extracted using the constant u 
method. We discard the results with no > 4 obtained from the g^ = extrapolation method. 



Constrained Fits to Aq and A3 

To gain more control over the extrapolations necessary to determine Ao{0), we combine them 
with a determination of ^3(0) and apply the constraint Ao{0) = ^3(0). We have previously 
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Figure 1 Comparison of linear and quadratic in 1/M extrapolations of ^o(O) ("(7^ = extrapo- 
lation" method). AQ{q^) has been fitted at each heavy mass value. The values of ^o(O) from the 
individual fits (squares) are then extrapolated linearly (solid curve) or quadratically (dashed curve) 
in 1/M to the B scale (diamond or cross). In the left hand plot, the extrapolated points have been 
displaced slightly for clarity. On the left a dipole form [no = 2 in equation (^)] has been used at 
each heavy mass, and on the right a quadrupole (ng = 4) form. 




Figure 2 AQ{(f') for the decay p'^l~vi obtained by the constant uj method. Squares denote 

points extrapolated from measurements around the charm scale in channels where oo is constant 
(or nearly constant). The curves shown are fits to equation ([l^ ) with values no = 1,5 (solid) and 
no = 2,3,4 (dotted) for the polar power. The three horizontal bars at = show the value 
and errors for Ao(0) when no = 2, and the burst point shows the value obtained from the = 
extrapolation method, also with no = 2. For both methods, the value of Ao(0) decreases as no 
increases. 
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Figure 3 Comparison of q'^ dependences of Aq and at (a) the charm (left) and (b) the bottom 
(right) scales. The squares are values for ^-nd the diamonds are values for A^. The charm plot 
in (a) uses our data with k = 0.129 for the heavy quark and with the light quarks extrapolated to 
'^crit = 0.14315(1) |2^]. The curves on the left hand plot are fits using the model of equations (16) 
and (ITtI) , with hq = 1 (solid) and no = 5 (dotted), as explained in the text. 



used this idea to perform constrained fits to and /° in 7r+/~z/; as well as Ti and 

T2 for B —>■ K*'~f. To fit both Aq and ^3 with "resonant" (pole or multipole) behaviour in 
g^, whilst maintaining the constraint Ao{0) = ^3(0) and agreement with the heavy quark 
symmetry requirements of equation (H) at g^ax' would require a model of the form, 

We find that such a model cannot describe our data. 

We were able to extract A3 (5'^) at the B scale for four momentum channels in our lattice 
simulation. The channel with the largest value of g^, 1 — (the zero recoil channel, — 
cannot be extracted directly for ^3), was too noisy to be extracted reliably: as we shall see, 
better quality lattice data allowing its extraction would be the easiest way to improve our 
results. Our results for A^^q"^) after extrapolation to the B scale in the constant uj method 
are included in table |l] (fourth column of table). 

We find that ^3 is positive for close to q^^.^^ in our measured data, with heavy quark 
masses around the charm scale. However, once extrapolated to the bottom mass, ^3(5^) is 
negative for close to g^ax- ^3 is a linear combination of Ai and A2, with positive and 
negative coefficients respectively, where the Ai contribution is most important at the charm 
scale, but Ai and A2 contribute more equally at the bottom scale. Moreover, heavy quark 
symmetry predicts that Ai decreases with increasing heavy quark mass, while A2 increases. 
In figure ^ we show the g^ dependence of Aq and ^3 at the two scales (the meaning of the 
curves in the charm-scale plot is explained later). 

The negative values of A^ found from our simulation prompt us to make two observations: 
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0.0 0.2 0.4 0.6 

m/M 

Figure 4 1/M extrapolation for A3 for ttie 0^1 channel, from the charm scale to the bottom 
scale. The squares are measured points for our four heavy mass values, which are extrapolated 
linearly (solid curve) or quadratically (dashed curve) in 1/M to the B scale (diamond or cross). 
The extrapolated points have been displaced slightly for clarity. 



1/M corrections to the heavy quark limit are large for ^3, as shown in figure ^. 
has the behaviour ^3 ~ MA2 ~ M^/^ in the limit M ^ 00, but in extrapolating from 
scales around the charm mass, we start from a region where the Ai contribution, with 
behaviour MAi ~ M^/^, is non-negligible. A3 exhibits the largest 1/M corrections we 
have found so far in several studies of processes involving heavy-to-heavy |32|, |33 



or 



heavy-to-light ||2^, ^ quark transitions. This is confirmed by using Ai and A2 values 



extracted previously by us |^ to calculate ^3 from the definition of equation (Q): 
these calculated values are given in the last column of table |l] as Ag^''^. We can also 
reproduce the coefficients found for the 1/M expansion of y430(ms/M)^''^ by using the 
definition of equation (^) to express them in terms of the 1/M expansion coefficients 
of Aie(mB/M)~i/2 g^^^ A2Q{mB/My/^ found in reference 



Resonant contributions alone, meaning dependences of the form given in equa- 
tions ([1^) and ([T5|) , which satisfy the constraint Ao{0) = ^3(0), are incompatible with 



the measured difference in sign of Aq and A3. The fact that around g^ax' ^0 is positive 
and A3 is negative, is a clear failure of pole dominance models to explain the be- 
haviour of the form factors in a region away from the physical pole. Sum rule P] - p3| 



and quark model [0, calculations, and the analysis of reference [|18|, also indicate 
that simple single pole fits cannot be made to all form factors. 

The negative values found for ^3 near q'^^^ are, in fact, still compatible with heavy quark 
symmetry constraints as M ^ 00. To make a combined fit for Aq and ^3, we use the 
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following dependences^]: 



Mq') = , (16) 

Mf) = ^o.(o)^if4^- (17) 



'1 — /m\ 



These forms satisfy the constraint Ao(0) = ^3(0) and are compatible with heavy quark 
symmetry provided that ^0,3(0) ~ M~'^°^^^'^ and c is independent of M. 

Because we had only four points for ^3(9^), covering a limited range of g^, we were unable 
to determine simultaneously the parameters ^9,3(0), c, mo and m^, in equations ([TBD and ([T7|). 
We therefore set = mo, as predicted by heavy quark symmetry in the infinite mass limit. 
In fact we believe that the l"*" B resonance should be heavier than the 0~ resonance and 
therefore that should be flatter in than Aq. By constraining both masses to be equal 
we make an error apparently at quadratic order in the expansion of the form factor A3. 
However, the fit parameter c can absorb quadratic differences and our errors in fact begin 
at 0{q^). 

This parameterisation fits our data well at the charm scale, as shown in figure |^a: the 
figure shows fitted curves for Aq and A3 obtained for tiq = 1 and uq = 5. Over a limited 
range of g^, changes in uq can be compensated by variations in the fitted mass parameter, 
mo. 

In figure we show the results of the combined fit at the B scale for several choices of the 
polar power no in the model of equations (|1^) and (|1^). The values obtained for ^40,3(0) are 



listed in table |] as one of the columns in the set labelled "constant u" . The table also gives 
^0.3(0) obtained from constrained fits for each of our heavy quark mass values, followed by 
an extrapolation of v4o,3(0) to the B scale (g^ = extrapolation method), including terms 
up to quadratic in the 1/M extrapolation. From the results of table |^ we can estimate the 
error due to the approximation m3 = mo. The largest error corresponds to the momentum 
channel — *• 1, which has the largest value, g^ = 17.5(2) GeV, used in the analysis for ^3. 
For this channel, taking m3 — mo ~ Aqcd ~ 0.5 GeV (in the D meson system the analogous 
mass difference is m£)^(i+) — m£)0(o-) = 2.423 GeV — 1.865 GeV = 0.558 GeV), the error is 7% 
for no = 1, ... , 5. 

In figure ^, the ^3 curves drop steeply as g^ approaches the pole at m| = mg. This looks 
like an artefact of using the same mass parameter for both Aq and ^3, when we expect m3 
to be heavier. Higher powers of no, while artificial, effectively increase the pole mass and 
make the g^ behaviour of A^ flatter away from the pole region. 

We find a range (0.16-0.35)lg for Ao(0). We do not include a 10% systematic error due 
to the chiral extrapolation and discretisation effects, as discussed in reference because 



this error is much smaller than both the statistical and systematic (owing to the different g - 
dependences assumed for Aq and ^3) errors noted above. Our range for Ao{0), (0.16-0.35)i^, 
is compatible with the value obtained in reference [0: Ao{0) = 0.24 ± 0.02. 

The constrained fits give the same result as fitting Aq alone above, but having data 
for A2, close to g^^^^ can help choose between different pole behaviours. Improved lattice 
simulations in the near future with smaller statistical and systematic errors will provide 
much better quality data for A3 close to g^ax- This will allow us to remove the restriction 
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These are motivated by Stech's form factor model |16 and by similar forms used recently by L Lei- 



louch p4] for the form factors /+ and in semileptonic B ^ n decay. 
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Figure 5 Combined fit to ^-nd ^3 at the B scale using the model of equations (16) and (17) 
for the solid and dotted curves respectively, with different values of the polar power uq as labelled 
on each plot. Squares are values for Ao(g^) and diamonds are values for ^3(^^)5 the same in all five 
plots. The three horizontal bars at = show the value and errors for ^0,3(0) = Ao{0) = ^3(0) 
in the combined fit (constant to method), and the burst point shows ^0,3(0) determined from the 
= extrapolation method. 
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^0,3(0) 




constant uj 

c mo/GeV x^/dof 



+7 
-6 


0.37: 


h6 
-5 


4.8l? 


0.31 


+8 
-6 


0.36: 


h6 
-5 


5.51^ 


0.34 


tl 


0.36: 


h6 
-5 


e.olg 


0.36 


+9 

-6 


0.36: 


h6 
-5 


6.611 


0.37 


+9 

-6 


0.36: 


h6 
-5 


7.11^ 


0.38 



Table 3 Values of the form factor Ao(0) = ^3(0), together with the fitted parameters c and ttiq, for 
p^l^vi from a combined fit of Aq and A3 using the model given in equations ( [I^ ) and (p!7|). 
Errors are statistical only. 



^2(0) 



Reference 



0.281^1^ 

0.24 ±0.56 

0.27 ±0.80 

0.38 ±0.18 ±0.04 

0.49 ±0.21 ±0.05 

0.721111^0 



this work 
APE "a" 
APE "b" m 



23 



m 



ELC "a" 



ELC "b" |2T 



22 



Table 4 Lattice determinations of ^2(0). The letters "a" and "b" denote results from two 
alternative methods of 1/M extrapolation (see [20| and [21 1 for details). 



= mo and will help to discriminate between the results for the different dependences 
reported in table ^and figure ^. For example, if A^i^q"^) does not decrease as — > g^ax' then 
lower values of Aq{0) will be preferred. 

In a previous paper ||3^, we showed that Ai(0) could be extracted using an assumed pole 



form for its dependence. Here we have extracted Aq{0) = ^3(0), which, by the definition 
of A^ in equation (^, makes possible a determination of ^2(0). This can help constrain 
fits to A2. Using a value Ai(0) = O.27I4II taken from our earlier results in reference ^0 



we findQ ^2(0) = 0.28lgl5, where in both cases the first errors are statistical and the second 
ones are systematic, due to discretisation effects and uncertainties in the chiral extrapolation 
together with the variation from changing no. Since ^2(gmax) > ^il^'max)' that A2 



has a more pronounced g dependence than Ai: this is consistent with sum rule |11] and 



quark model JTB], [T^ calculations, and with the analysis presented in fT^. Our result for 



^42(0) is the most precise obtained from the lattice to date: it is compared with previous 
lattice determinations in table ^ 

Unfortunately, the quality of the numerical results for Aq{0), obtained here, and for T2 
in our previous work [^, does not allow us to calculate the ratio 2T2^^^*(0)/yl^~*''(0) of 
hadronic form factors in equation (|^). In reference [^, we found two sets of values for ^2(0), 
0.15l6 and 0.261^ depending on the g^ behaviour (pole or constant, respectively) assumed 
for the form factor T2(g^). Some authors ^j, have suggested that the ratio / is close to 
one. This, together with the results for Aq presented here, would support low values for T2(0) 

^We quote A2{0) obtained from ^3(0) with no — 2. The central value of ^2(0) ranges from 0.24 for 
no = 1 to 0.30 for uq = 5. This range of values is smaller than the size of the statistical errors quoted above. 
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and therefore a pole- type behaviour for its dependence. In reference |3^, by studying the 
ratio Ai/2iT2 (predicted to be 1 in the heavy-quark hmit) and the form factor Ai, we also 
found arguments supporting a pole-type behaviour for T2. 

This study, combined with our previous studies of the form-factors Ai, A2 and V in 
5° ^ P+/-Z// decays 0], Tl and T2 in i? K*-^ decays and /+ and f in 5° ^ 7r+/"z/i 



decays , will enable us in the immediate future to perform a global analysis of the heavy- 
to-light decays tt^I^iji, B —>■ K*'~f and B^ p~^l~ui. Such a study will serve to 
check the validity of different theoretical models, for example the one recently suggested by 
In conjunction with theoretical input, such as the constraints on and /° 

m 



Stech 
B^ - 



n+t 



m 



ui derived from a dispersion relation analysis in 
calculations reduce the systematic uncertainties related to the q 
factors. 



it could also help lattice 
^ dependence of the form 



Conclusions 

We have determined a range of values for Aq{0) from the lattice for the first time. The 
range we find is: (0.16-0.35)tg. As data on both Aq and A3 improves, it could become an 
important ingredient in the determination of \Vub\- 

From fits to Aq alone, it is hard to dismiss high values for Ao{0) arising from an assumed 
single pole behaviour of the form factor. Information from ^3, although not of high quality, 
can help determine Ao{0) by disfavouring such high values. With our present data we cannot 
rule out the higher values, but with better quality data, in particular A^i^q"^) values with 
closer to ^max) ^ constrained fit to Aq and A3 could give more reliable answers (as was found 
for the form factors and for B^ tt^I^ui in |^). This would allow us to reduce the 
range of values quoted for Aq{0). 

A3 itself, for q^ near g^^x; a large dependence on the heavy quark mass, changing 
from positive to negative as one scales from the charm mass to the bottom mass. Because 
A^i^q'^) at the B scale is negative when Ao(g^) is positive at the same g^, close to g^ax; 
cannot be fitted solely by resonant terms. Hence, our lattice measurements bear out the 
predictions of quark models. 

Combining the results for Aq{0) with earlier results for Ai allows us to determine ^2(0) = 
0.28lg^5. Combining this with improved lattice measurements of v42(g^) for g^ near q'^^^^ 
could allow the determination of the full g^ dependence of A2. 

Future lattice simulations, with improved statistical and systematic errors, would allow us 
to take m3 7^ mo in the fits to equations (|TB|) and (P^). We believe that in these circumstances 
the most physically motivated polar power, no = 1, would be preferred. We expect that the 
results obtained would coincide with those presented here, where we have set = mo (as 
predicted by heavy quark symmetry) and have limited ourselves to a g^-region away from the 
pole where higher polar powers effectively increase the pole mass and make the g^ behaviour 
of A3 flatter. 
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